void AEPD_Fit(double* x, int N, double* params);
void fast_AEPD_Fit(double * x, int N, double * params, struct DeviationStruct *D);
double AEPD_PDF(double* x, int N, double theta, double sigma, double kappa, double alpha);
double AEPD_CDF(double x,double theta, double sigma, double kappa, double alpha);
double AEPD_Inv(double Percentile, double theta, double sigma, double kappa, double alpha);
void AEPD_Stat(double theta, double sigma,  double kappa, double alpha, double* ret);
//void AEPD_Stat(double theta,double sigma,  double kappa, double alpha, double ret[2]);


